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The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has performed systematic 
measurements of <j> meson production in the K + K~ decay channel at midrapidity in p + p, d+Au, 
Cu+Cu and Au+Au collisions at y/s NN —200 GeV. Results are presented on the <f> invariant yield 
and the nuclear modification factor Raa for Au+Au and Cu+Cu, and RdA for d+Au collisions, 
studied as a function of transverse momentum (1 < pr < 7 GeV/c) and centrality. In central and 
mid-central Au+Au collisions, the Raa of cj> exhibits a suppression relative to expectations from 
binary scaled p+p results. The amount of suppression is smaller than that of the n° and the r\ 
in the intermediate pr range (2-5 GeV/c), whereas at higher pr the <fr, ty° and n show similar 
suppression. The baryon (protons and anti-protons) excess observed in central Au+Au collisions at 
intermediate pr is not observed for the <j> meson despite the similar mass of the proton and the <f>. 
This suggests that the excess is linked to the number of constituent quarks rather than the hadron 
mass. The difference gradually disappears with decreasing centrality and for peripheral collisions 
the Raa values for both particles are consistent with binary scaling. Cu+Cu collisions show the 
same yield and suppression as Au+Au collisions for the same number of iVp ar t. The RdA of <j> shows 
no evidence for cold nuclear effects within uncertainties. 

PACS numbers: 21.65. Jk,25.75.Dw 



I. INTRODUCTION 



Measurements of hadron spectra from p + p and 
nucleus-nucleus collisions at RHIC provide a means to 
study the mechanisms of particle production and the 
properties of the medium formed in relativistic heavy ion 
collisions. At low transverse momentum, pt < 2 GeV/c, 
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where the bulk of particles are produced, hadron produc- 
tion is governed by soft processes characterized by low 
momentum transfer. The particle yields and the evolu- 
tion of the interacting system are successfully described 
within the framework of thermal and hydrodynamical 
models [lH|. 

At high transverse momentum, px > 5 GeV/c, hard 
scattering processes become the dominant contribution. 
Due to the large momentum transfer involved, the 
parton-parton scattering cross sections are amenable to 
perturbative QCD (pQCD) description and hadron pro- 
duction can be calculated using initial state parton dis- 
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tribution functions and final state fragmentation func- 
tions. Modifications to the hadron yields are expected 
in nucleus- nucleus collisions due to the interaction of 
the scattered parton with the hot and dense medium 
formed In the absence of interaction with the 

medium the hard scatterings and the resulting hadron 
yields should scale with the number of binary nucleon- 
nucleon collisions (N co ii), whereas in the medium the 
yields are suppressed ("jet-quenching" @) due to par- 
ton energy loss through gluon bremsstrahhmg. High- 
er hadron suppression consistent with this scenario has 
been discovered in Au+Au collisions at RHIC [ltJ-Qjl]. 
The same suppression by a factor of ~ 5 is observed for 
7T° and r) production whereas direct photons which do 
not interact with the medium, follow the expected bi- 
nary scaling (l3| . Single electrons originating from the 
semi-leptonic decays of mesons containing heavy quarks 
(charm and bottom) exhibit a large suppression at high 
Pt, similar within the experimental uncertainties to that 
of 7T° and ij, presenting a challenge for the bremsstrahlung 
explanation flij ]. 

At intermediate transverse momentum 2 < 
Pt (GeV/c)< 5, suppression of binary scaled pro- 
duction is observed for light ir° and r\ mesons but not 
for protons and a nti-p rotons in mid-central and central 
Au+Au collisions [l5[ . The p/ir and p/ir ratios increase 
with centrality and exceed the values measured in p + p 
by a factor of 3-5 in the most central collisions. A dif- 
ferent suppression pattern between baryons and mesons 
is also observed for strange hadrons, A and Kg [l6l 
These baryon/meson differences in suppression are 
inconsistent with the picture of hadron production 
through hard-scattering followed by partonic energy loss 
in medium and hadronization in vacuum according to 
the universal fragmentation functions. This poses the 
question whether hard-scattering is the dominant source 
of baryon production at intermediate pt- Studies of jet- 
like dihadron correlations in Au+Au collisions [l8|, Tl9j 
imply nearly equal importance of the jet fragmentation 
as a production mechanism for mesons and baryons, 
except for the most central collisions. The interpretation 
of the baryon non-suppression results requires therefore 
another particle production mechanism in addition to 
jet fragmentation at intermediate pt- 

There have been attempts to describe the different be- 
havior of baryons and mesons through the strong radial 
flow that boosts particles with larger mass to higher 
Pt [20l , I2H or through the recombination of soft and 
hard massive partons [22T424T ] . or through the interplay 
of the jet quenched hard component and phenomenolog- 
ical soft to moderate pt baryon junction component [25j 
or through the QCD color transparency of higher-twist 
contributions to inclusive hadroproduction cross sections, 
where baryons are produced directly in short distance 
subprocess [261 ]. Although several models reproduce pt 
spectra, particle ratios and elliptic flow for different 
hadrons reasonably well, the relative contributions from 
the different processes are difficult to infer. 



The (ft meson is a very rich probe of the medium formed 
in heavy ion collisions, because it is sensitive to several 
aspects of the collision, including strangeness enhance- 
ment and chiral symmetry restoration, as well as energy 
loss and the nuclear modification factor [2TM3ll j , which is 
the focus of this paper. Due to its small inelastic cross 
section for interaction with non-strange hadrons [27|, [Hj] , 
the 4> meson is less affected by late hadronic rescatter- 
ing and reflects better the initial evolution of the system. 
Being a meson with a mass comparable to that of the pro- 
ton, it is interesting to see how the <f> meson fits within 
the meson/baryon pattern described previously; being a 
pure ss state, it puts additional constraints on the energy 
loss and recombination models. 



This paper presents systematic PHENIX measure- 
ments of 4> meson production via the K + K~ decay chan- 
nel at v /s^7 =200 GeV, including first PHENIX results 
in p+p, e?+Au and Cu+Cu collisions and new results in 
Au+Au collisions. The latter have much higher statis- 
tics and a finer centrality binning in c omp arison to the 
previously published PHENIX results (Hf. The results 
benefit from three different techniques involving differ- 
ent levels of kaon identification in the analyses. These, 
combined with the high statistics of the analyzed data 
samples, allow for the extension of the pt range of the 
measurements up to pt — 7.0 GeV / c in all collision sys- 
tems. The higher pt reach and the higher precision of 
the data allow for sharper conclusions with respect to 
earlier results [28|, |3(J. The Cu+Cu measurements arc 
complementary to those on Au+Au and allow the study 
of nuclear effects with different nuclear overlap geometry 
for the same N vax t and with smaller N part uncertainties 
for N paTt < 100. 



The measurement of the <fi meson production in d+Au 
collisions is important for understanding cold nuclear 
matter effects which are of interest by themselves and 
are also essential for the interpretation of heavy ion col- 
lisions. As shown in [33j . in the intermediate pt range, 
charged pions are practically not enhanced in compari- 
son to the binary scaled p + p yield, whereas protons and 
anti-protons exhibit some enhancement of ~ 30% in the 
most central collisions. The mechanism of multiple soft 
re-scattering of partons in the initial state which is usu- 
ally invoked as the origin of the Cronin effect does not 
explain this meson/baryon difference. One possible ex- 
planation comes from recombination models |34j in which 
baryons gain higher transverse momentum from recom- 
bination of three quarks in the final state in comparison 
to mesons consisting of only two quarks. Measurement of 
the Cronin effect for the <fi mesons can provide additional 
constraints for the models that try to explain these cold 
nuclear effects. 
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II. EXPERIMENTAL SET-UP AND DATA 
ANALYSIS 

We report on the measurements of <f> mesons at midra- 
pidity in the K + K~ decay channel in p+p, d+Au, Cu+Cu 
and Au+Au collisions at y/s NN —200 GeV using data 
collected by the PHENIX experiment during the 2004, 
2005 and 2008 physics runs. A detailed description of 
the PHENIX detector can be found elsewhere [35]. The 
measurements were done using the two PHENIX cen- 
tral arms each covering 90° in azimuth at midrapidity 
(| ry | < 0.35). The tracking of charged particles and the 
measurement of their momentum with typical resolution 
of Sp/p = 0.7© l.l%p [GeV/c] are performed using the 
drift chambers and the first layer of the pad chambers 
(PC). In order to reduce the background at high px, 
tracks are required to have a matching confirmation in 
the third layer of PC or the electromagnetic calorime- 
ter (EMCal). Kaons are identified using the time-of- 
fiight detector (TOF), which covers approximately 1/3 
of the acceptance in one of the central arms. With a 
time resolution of ~ 115 ps, the TOF allows for clear 
7r/K separation in the range of transverse momentum 
from 0.3 GeV/c to 2.2 GeV/c using a 2a p^-dependent 
mass-squared selection cut as described in [28[ . 

The beam-beam counters (BBC) and zero degree 
calorimeters (ZDC) are dedicated subsystems that de- 
termine the collision vertex along the beam axis (z vtx ) 
and the event centrality, and also provide the minimum 
bias interaction trigger. Events are categorized into cen- 
trality classes using two-dimensional cuts in the space of 
BBC charge versus ZDC energy [36[ for Au+Au colli- 
sions or only by the amount of charge deposited in the 
BBC 0, |37l for d+Au and Cu+Cu collisions. 

In any particular event one cannot distinguish between 
kaons from <f> decays and other kaons, so the <p meson 
yields are measured on a statistical basis. In each event, 
all tracks of opposite charge which pass the selection cri- 
teria are paired to form the invariant mass distribution. 
This distribution contains both the signal (S) and an in- 
herent combinatorial background (B). To maximize the 
statistical significance and the px reach of the measure- 
ments we use three different track selection techniques: 
"no PID" , in which all tracks are assigned the kaon mass, 
but no TOF information is used, and "one kaon PID" or 
"two kaons PID" , in which one or both tracks are iden- 
tified as kaons in the TOF. 

Table U lists for each collision system and for each anal- 
ysis technique the number of analyzed minimum bias 
events in the vertex range |z v tx| < 30 cm, the accessi- 
ble pt range and the range of the S/B ratio. 

The raw yields of the <fi are obtained by integrating the 
invariant mass distributions in a window of ±9 MeV/c 2 
around the ci mass after subtracting the combinatorial 
background. In the analysis of Au+Au, Cu+Cu and 
d+Au data, the combinatorial background is estimated 
using an event-mixing technique. The details of the 
method are given elsewhere [28j. In the "no PID" anal- 



TABLE I: Collision species, number of analyzed minimum 
bias events, accessible pr range and typical range of S/B ratio 
for the different (f> — > K + K~ analyses. 



Species N [10 9 ] p T [GeV/c 



S/B 



Technique 



P+P 
d+Au 
Cu+Cu 

Au+Au 



1.50 
1.44 


0.9- 
1.3- 


-4.5 
-7.0 


1/9 
1/76 


-1/2 
- 1/3 


"one kaon PID" 
"no PID" 


1.69 


1.1- 


-7.0 


1/245 


- 1/12 


"no PID" 


0.77 
0.78 


1.1- 
1.9- 


2.95 
-7.0 


1/91 
1/205 


- 1/9 

- 1/24 


"one kaon PID" 
"no PID" 


0.72 
0.82 


1.1- 
2.45 


3.95 
-7.0 


1/19 
1/385 


- 1/2 

- 1/32 


"two kaons PID' 
"no PID" 
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FIG. 1: (Color online) Invariant mass distributions obtained 
with the "two kaons PID" and "no PID" methods in Au+Au 
collisions after subtraction of the combinatorial background 
estimated using the event-mixing technique. Plot on the top 
corresponds to integrated pr range whereas plot on the bot- 
tom is for the range 2 < pr (GeV/c)< 3. The "no PID" spec- 
trum is fit to the sum of a Breit-Wigner function convolved 
with a Gaussian function to account for the </> signal, and a 
polynomial function to account for the residual background. 



ysis, a significant residual background remains in the 
subtracted mass spectra because the mixed-event tech- 
nique does not account for the abundant correlated pairs 
from other particle decays (K° — > 7r + 7r~, A — > pn~ , 
p — > 7r + 7r~, u) — > ir°Tr + ir~ etc.). In the "one kaon 
PID" analysis the residual background is considerably 
smaller [38[ while in the "two kaon PID" method the 
background is negligible. Examples of subtracted mass 
spectra obtained in Au+Au collisions with the "two kaon 
PID" and "no PID" techniques arc shown in Fig. [TJ The 
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S/B ratio depends on the collision system, the analysis 
technique, the <j) transverse momentum and the central- 
ity. The typical ranges of the S/B values for each colli- 
sion system and each analysis technique are summarized 
in Table |U 

The total combinatorial background in p+p |38{ . as 
well as the residual background in d+Au, Cu+Cu and 
Au+Au analyses were estimated by fitting the mass spec- 
tra with the sum of a Breit-Wigner mass distribution 
function convolved with a Gaussian experimental mass 
resolution function to account for the (f> signal, and a 
polynomial function to account for the background. The 
typical experimental mass resolution for the (f> meson was 
estimated to be ~ 1 MeV/c 2 using Monte-Carlo studies 
based on the known momentum resolution of the tracking 
system and time resolution of the TOF. To describe the 
background a second order polynomial was used in most 
analyses, except for the Au+Au "no PID" case where 
a third order polynomial was used. Figure [1] shows an 
example of the fits. 

The <f> meson invariant yield in a given centrality and 
Pt bin is obtained by: 
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o 
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Au+Au, Min. bias ■ two kaons pid 
• no PID 
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FIG. 2: (Color online) Top: Ratios of yields obtained with 
"no PID" and "one kaon PID" (in p+p collisions), or "no 
PID" and "two kaons PID" (in Au+Au collisions) to fits to 
the combined spectra. Bottom: comparison of the acceptance 
and reconstruction efficiency for the three different analysis 
techniques. 



1 d 2 N A^Cbias (1) 

2np T dprdy 2TTpTN cvt e vcc e cm bcdBKK^PT^y ' 

where iV cv t is the number of analyzed events in the cen- 
trality bin under consideration, e rcc corrects for the lim- 
ited acceptance of the detector and for the <j> meson re- 
construction efficiency, e e mbed accounts for the losses in 
reconstruction efficiency due to detector occupancy in 
heavy ion collisions, £>kk is the branching ratio for <j> — ¥ 
K + K _ in vacuum, is the raw <fi yield measured in 
the given bin, C bias = ^ftf* , where e^| c and e^ BC 
are the beam-beam trigger efficiencies for minimum bias 
and <f> events res pec tively. This Cbi as correction is equal 
to 0.69 for p+p [39j and varies from 0.92 to 0.85 as we 
go from peripheral to central d+Au collisions [4(j. In 
Au+Au and Cu+Cu collisions the minimum bias trigger 
is inefficient only for very peripheral collisions (central- 
ity > 92.2% for Au+Au and > 94% for Cu+Cu). For 
all other centralities, - 92.2% (0 - 94%) for Au+Au 
(Cu+Cu), there is no trigger bias and Cbias = 1. In p+p 
the invariant differential cross section at midrapidity is 
related to invariant yield as = a™ el x £^ y , 

where a^E , =42.2 mb. 

The detector occupancy related loss (1 — e em bed) is cal- 
culated by embedding simulated K + K~ pairs into real 
events. It varies from 1% in peripheral to 29% (7%) in 
the most central Au+Au (Cu+Cu) collisions without a 
significant px dependence. 

The acceptance and reconstruction efficiency (e roc ) of 
the 4> meson, determined using a full GEANT simulation 
of the PHENIX detector, is shown in the bottom panel 
of Fig. [5] for different analysis techniques. There are very 
large differences, reaching more than one order of mag- 
nitude between the three cases. In spite of that, the 



invariant yield spectra obtained from the different tech- 
niques are in good agreement as demonstrated in the top 
panel of Fig. [21 which shows the ratios of yields obtained 
with "no PID" or with "one kaon PID" ("no PID" or 
"two kaons PID") techniques in p+p (Au+Au) to a fit 
performed to the combined data sets. This agreement 
implies good control over the systematic uncertainties 
which are quite different in the three cases and provides 
confidence on the robustness of the experimental results. 

The results from measurements at low px using "two 
kaons PID" (in Au+Au collisions) and "one kaon PID" 
(in p+p and Cu+Cu) are combined with the independent 
"no PID" measurements at intermediate and high pt to 
form the final px spectra. The measurement in d+Au 
is performed using the "no PID" technique only. The 
invariant mass spectra obtained with "one kaon PID" 
or "two kaon PID" methods are subsamples of the "no 
PID" distribution. Therefore results obtained with dif- 
ferent methods can not be directly averaged. In the final 
spectra the transition between different techniques occurs 
at p T = 1.3 GeV/c in p+p, p T = 2.2 GeV/c in Au+Au 
and at px — 3.2 GeV/c in Cu+Cu collisions in order to 
obtain the smallest combined statistical and systematical 
uncertainties for the points. 

Systematic uncertainties on the <f> invariant yield are 
grouped into three categories: type A (point-to-point un- 
corrected), which can move each point independently; 
type B (point-to-point py-correlated), which can move 
points coherently, but not necessarily by the same rela- 
tive amount; type C (global), which move all points by 
the same relative amount. The main contribution to the 
systematic errors of type A is the uncertainty in the raw 
yield extraction N<j, of 6-25%. Error of type B is dom- 
inated by uncertainties in reconstruction efficiency e roc 
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FIG. 3: (Color online) Invariant pr spectra of the <f) meson for different centrality bins in Au+Au, Cu+Cu, ei+Au and p-\-p 
collisions at y/s NN =200 GeV. The statistical and systematic uncertainties are smaller than the size of the symbols. The spectra 
are fitted to exponential and Tsallis [4ll-|43T] functions shown by the dashed and solid lines, respectively. 



of 5-9%, embedding corrections e em bed of 1-7% and mo- 
mentum scale of 1-5%. The main contributions to the 
type C errors are the uncertainties in normalization for 
the p+p (d+Au) cross section equal to 9.7% (7.8%) and 
in branching ratio Bkk of 1.2%. 



III. RESULTS AND DISCUSSION 

Figure [3] shows the fully corrected <j) invariant yield 
as a function of pr measured in p+p, d+Au, Cu+Cu 
and Au+Au collisions at y/s^" =200 GeV. The spectra 
are scaled by arbitrary factors for clarity and fitted to 
exponential and Tsallis [4l| - |43j functions shown by the 
dashed and solid lines, respectively. We used the Tsallis 
function adapted to the form (38|: 

1 d 2 N _ 1 dN (n-l)(n-2) 
2ir dydpr 27r dy (nT + m^(n — l))(nT + m^) 

\ nT + m J 



where n and T are free parameters, tut = 

\Jpr + m ^> an d m is the mass of the 4> meson. The 
spectral shapes for all collision systems and centralities 
are well described by the Tsallis function, while the ex- 
ponential fits underestimate the <f> meson yields at high 
Pt where the spectra begin to exhibit the power law be- 
havior expected for particles produced in hard scattering 
processes. For p+p collisions the departure from expo- 
nential shape occurs at ~4 Gev/c. For all centralities 
in Au+Au collisions the departure occurs at somewhat 
larger px, which suggests a larger contribution of soft 
processes to the <f> meson production up to 4-5 GeV/c. 
Such behavior of the spectral shap es is in agreement with 
recombination models [22hl23 . l44 - |4r| predicting pt spec- 
tra for different hadronic species based on the number 
and flavor of constituent quarks. At low transverse mo- 
mentum, we do not observe a large change in the slopes of 
the spectra from central to peripheral collisions, support- 
ing the expectation for smaller radial flow in mesons 
compared to other hadrons. 

The large pt reach of the results presented here allows 
for the study of medium-induced effects on (j) meson pro- 
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duction at intermediate and high px using the nuclear 
modification factor: 



Rab(pt) = dN AB /((N co n) x dN pp ), 



(3) 



where dN ab (dN pp ) is the differential <j> yield in nucleus- 
nucleus (p+p) collisions and (N co n) is the average number 
of nuclear collisions in the centrality bin under consider- 
ation [ll], . The latter is determined solely by the 
density distribution of the nucleons in the nuclei A and 
B and by the impact parameter and is calculated using 
the Glauber formalism j47| . Deviation of Rab from unity 
quantifies the degree of departure of the A+B yields from 
a superposition of incoherent nucleon-nucleon collisions. 
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FIG. 4: (Color online) Top: Raa vs. pr for 0, 7r°, 77, 
(K + +K+) and (p + p) in central Au+Au collisions. Mid- 
dle: Raa vs. pr for <f> and it in 10-20% mid-central Au+Au 
collisions. Bottom: -Raa vs. px for (f>, and p + p in 60-92% 
and for tt° in 80-92% peripheral Au+Au collisions. Values 
for (K++K-), (p + p), tv° and 77 are from [H, [M • The 
uncertainty in the determination of {N co n) is shown as a box 
on the left. The global uncertainty of ~ 10% related to the 
p+p reference normalization is not shown. 



Figure @] shows a comparison of the Raa for 
7T° from Ref . 14911 . proton and kaon from Ref. [33 
77 from Ref. [50], all measured in Au+Au collisions at 
^s NN —200 GeV. The 4> meson exhibits a different sup- 
pression pattern than that of lighter non-strange mesons 



and baryons. For central collisions (top panel) the 4>'s 
Raa shows less suppression than 7r° and 77 in the inter- 
mediate pr range of 2 < pr (GeV/c)< 5. At higher 
Pt values, pt > 5 GeV/c, the (/>'s Raa approaches and 
becomes comparable to the 7r° and 77 Raa- These two 
features remain true for all centralities up to the most 
peripheral collisions as displayed in the bottom panel of 
Fig. 2] (see also Fig. [5]). The panel shows that the ir° is 
slightly suppressed (at the level of ~ 20%) in peripheral 
Au+Au collisions whereas the <j> is not suppressed. The 
kaon data cover only a very limited range at low pt but 
in this range they seem to follow the Raa trend of the (f> 
better than that of the 7r° and 77 for central Au+Au colli- 
sions. The comparison with baryons, represented in Fig- 
ure 2] by the protons and anti-protons, shows a different 
pattern. For central collisions, the protons show no sup- 
pression but rather an enhancement at px > 1.5 GeV/c, 
whereas the <f> mesons are suppressed. This difference be- 
tween <j> mesons and protons gradually disappears with 
decreasing centrality and for the most peripheral colli- 
sions the Raa of <p and (anti)protons are very similar as 
demonstrated in the bottom panel. 

The results presented here are in agreement with the 
previous PHENIX results (28[, which were based on the 
2002 RHIC run, within the relatively larger uncertain- 
ties of the latter. The use of different analysis techniques 
and the larger Au+Au data sample of the 2004 run re- 
sulted in a higher precision and a larger px reach of Raa 
that allowed to unveil the different behavior of the <j) me- 
son, i.e. less suppression than ir° but more suppression 
than baryons, in the intermediate px range. Our results 
differ from the ones recently published by the STAR Col- 
laboration [29|, [33| which show that in Au+Au collisions 
Raa is consistent with binary scaling in the intermediate 
px region whereas Rcp shows considerable suppression. 
This difference is traced down to the almost factor of two 
higher invariant px yield in the STAR experiment [29|, [3(| 
in Au+Au collisions, compared to our results presented 
in Fig. [31 whereas in p + p both experiments are in rea- 
sonably good agreement. 

Figure[5]compares the Raa of <j> in Au+Au and Cu+Cu 
in two centrality bins which correspond approximately 
to the same number of participants in the two sys- 
tems. Figure [S] shows the Raa of the <j> integrated for 
px > 2.2 GeV/c in Cu+Cu and Au+Au collisions versus 
Npaxt. Under these conditions, there is no difference in 
the Raa of (f> between the two systems indicating that 
the level of the suppression, when averaged over the az- 
imuthal angle, scales with the average size of the nuclear 
overlap, regardless of the details of its shape. This behav- 
ior has been observed in other measurements, such as the 
Raa of the ir°. The 7r° suppression data in Au+Au and 
Cu+Cu taken from Ref. [l2j)|49| are also shown in Fig. [S] 
for comparison. The similarity of the Raa of <f> in the two 
colliding systems implies that the features discussed pre- 
viously for Au+Au in the context of Fig. namely that 
the <j) exhibits an intermediate suppression between pions 
and baryons, remain valid also in the Cu+Cu system. 
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FIG. 5: (Color online) Top: Raa vs. pr for </> and 7r° for 
30-40% centrality Au+Au and 0-10% centrality Cu+Cu col- 
lisions. Bottom: Raa vs. pr for cf> and n° for 40-50% central- 
ity Au+Au and 10-20% centrality Cu+Cu collisions. Values 
for 7T° are from [13, IS]. The uncertainty in the determination 
of (JVcoii) is shown as a box on the left. The global uncertainty 
of ~ 10% related to the p+p reference normalization is not 
shown. 



Our data disfavor radial flow as the dominant source 
for the particle species dependence of the suppression fac- 
tors at intermediate p?, since the proton and <f> Raa fac- 
tors differ by a factor of ~ 2, in spite of their similar mass 
(m p ~ nifj,), whereas the kaon and <j> show similar Raa 
factors although their masses differ by almost a factor of 
two (to^ ~ 2771k)- 

Recombination models [22T - I241 |43 - |46| qualitatively ex- 
plain the larger yield of baryons compared to mesons at 
intermediate pt by the higher gain in pt which comes 
from recombination of three quarks for baryons rather 
than two quarks for mesons. The same framework can 
be used to interpret the difference in suppression fac- 
tors for 7T° and 4> mesons. For 7r° production in the 
Hwa and Yang model (46} the contribution from the re- 
combination of thermal (T) and shower (S) partons be- 
comes comparable to that of the recombination of TT 
partons already at pr ~ 3 GeV/c. For the <j> however, the 
strangeness enhancement feeds preferentially the thermal 
quarks. Soft processes dominate over hard processes in 
a wider pt range and consequently the TT component 
remains dominant up to px ~ 6 GeV/c for the <t> produc- 
tion [45j ■ The Raa of <j> becomes similar to that for ir° at 
Pt > 5 — 6 GeV/c where the contribution from fragmen- 
tation partons becomes significant for both particles. It 
is interesting to note that the r\ follows closely the 7r° in 
spite of its sizable (~ 50%) strangeness content [5l[ . 
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Cold nuclear matter effects can also contribute to 
the differences in hadron suppression factors in A+A 
collisions. Figure [7] compares the RdA for <j> and 7r° 
from Ref. 48|, and protons from Ref. [33[ for central 
(top panel) and peripheral (bottom panel) d+Au colli- 
sions. For both centralities, the RdA for (j) and tt° are 
similar indicating that cold nuclear effects are not re- 
sponsible for the differences between <f> and tt° seen in 
Au+Au and Cu+Cu collisions. The proton's RdA ex- 
hibits an enhancement for px=2-A GeV/c, usually asso- 
ciated with the Cronin effect (52T457T ] . whereas the RdA 
for <fi indicates little or no enhancement. The lack of 
Cronin enhancement is also seen in the n° data (48j 
shown in Fig. [7] and has also been observed for other me- 
sons in central and mid-central d+Au collisions at ^s NN 

=200 GeV m, m, HI] ■ 

IV. SUMMARY AND CONCLUSIONS 

We have measured <f> meson production at midrapidity 
via the K + K~ decay channel in p+p, d+Au, Cu+Cu and 
Au+Au collisions at ^s NN =200 GeV. Invariant px spec- 
tra and nuclear modification factors have been presented 
over the px range of 1 < px < 7 GeV/c for different 
centralities. 

The <j) meson exhibits a different suppression pattern 
compared to lighter mesons (it and 77) and baryons (pro- 
tons and anti-protons) in heavy ion collisions. For all 
centralities, the (j) meson is less suppressed than tt° and 
•q in the intermediate pr range (2-5 GeV/c) whereas at 
higher px, 4>, 7r° and r\ show similar suppression values. 
The available kaon Raa data seem to follow the Raa 
trend of the <fi. The comparison with baryons shows that 
in central Au+Au collisions the latter are enhanced with 
respect to binary scaling whereas the 4> meson is sup- 
pressed, but this difference gradually disappears with 
decreasing centrality and for peripheral collisions the 
baryons and the (f> meson have very similar Raa values 
consistent with binary scaling. 

The same features are observed in Cu+Cu collisions 
between the 4> and it . The 4> meson invariant px spectra 
in Au+Au and Cu+Cu collisions for similar -/V part values 
exhibit similar shape and yield over the entire px range 
of the measurement within the statistical and systematic 
uncertainties. This indicates that the production and 
suppression of the <f> meson, when averaged over the az- 
imuthal angle, scales with the average size of the nuclear 
overlap region, regardless of the details of its shape. 

Cold nuclear effects cannot account for the observed 
differences. For all centralities, the <fi's RdA in d+Au col- 
lisions is consistent with binary scaling in agreement with 
other mesons. No meson species dependence is observed 
in RdA within uncertainties. 



The observed features at intermediate px in Au+Au 
and Cu+Cu collisions are qualitatively consistent with 
quark recombination models [22T4241 I3j46j . which are 
also supported by <p elliptic flow measurements [29|, [3l[ . 
The systematic set of measurements presented here pro- 
vides further constraints to these models. The similar- 
ity between the suppression patterns of different mesons 
at high pt favors the production of these mesons via jet 
fragmentation outside the hot and dense medium created 
in the collision. Complementary jet correlation measure- 
ments involving (f> mesons as a trigger as well as exten- 
sion of the kaon data to higher px would be desirable 
to provide further insight into the <f> meson production 
mechanism. 
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